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Abstract 

Dissipative Lagrangians and Hamiltonians having Coulomb, viscous and quadratic damping, 
together with gravitational and elastic terms are presented for a formalism that preserves the 
Hamiltonian as a constant of the motion. Their derivations are also shown. The resulting 
L's and H 's may prove useful in exploring new types of damped quantum systems. 
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Of the Lagrange-Hamiltonian formalisms that incorporate frictional forces, the formulation 
by Schuch pQ is closely analagous to a conservative system. This approach has had a number of 
advantages, especially for describing quantum systems. As expected, the Lagrangian L produces 
the correct equation of motion. The Hamiltonian H is a constant of the motion. H also has 
a physical interpretation of being the total energy for the system. ' The problem of shrinking 
Heisenberg uncertainty, often seen with canonical attempts at dissipative quantization (UElElllj, 
can be resolved. In addition, the range of candidates for H is greatly narrowed from an infinite 
set of possibilities, e.g. p3 E]. For convenience, the formalism developed by Schuch will be 
referred to as H exp , where exp represents an expanding coordinate space. 

The existence of H exp does not necessarily rule out the usefulness of alternate forms. Schuch 
Ej has also shown that well known formalisms that do not preserve H as a constant of the 
motion such as Caldirola [Jj and Kanai 0, can be converted to and from fl exp via canonical 
transformation. There is a similar relationship for the non-linear Schrodinger equation, (NLSE), 
EI- These other formalisms and H exp are closely related because they share similar non- 
canonical transformations that connect the physical and transformed coordinate spaces. 

Given its promising characteristics, new instances are presented to show that the H exp for- 
malism can be broadened to include additional combinations of friction and potential. Some of 
the more important details concerning H exp are discussed briefly. 

To obtain L's and H's that have the desired properties, it is necessary to determine relations 
that connect the physical space coordinate q, to a transformed Q space [TJ, 

Q = F(q,t,£o,ti, ■ ■ ■ ,4>0,<f>i, ■ ■ ■) (!) 

where t is time and & and (pi are the coefficients of friction and potential. Note that within 
Q-space explicit time dependence is absent. 

* E-mail address: cesmith@vt.edu 

t Depending on the initial conditions that are specified. 
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The resulting L's and LTs for a prospective system with frictional forces must also satisfy a 
set of formalisic rules which are listed below. 

1. All L(q,q,t)'s, L(Q,Q)'$, H(p,q,t)'s, and H(P,Q)'s have the units of energy. 

2. The L for the system has an invariance property whereby both SL(q, q, t) = and SL(Q, Q) 
= correctly produce the physical space equation of motion. * 



3. The physical and transformed momenta p and P are defined as p — mq and P = mQ 

dL(Q,Q) 
dQ 



respectively. In the absence of expressed dissipative forces, as in Q-space, P 

dq 



also holds. If friction is expressed explicitly, as in real space, p ^ 8 -L(<;,q,t) m general! 



4. H(P, Q) is found via the Lengendre transform H = PQ — L where neither L = L(Q, Q) 
nor H(P, Q) contain t explicitly. If there is explicit time dependence in the connection 
between Q and q, H(p, q, t) ^ pq — L(q, q,t)\ = 

5. The H (p, q, t) that is obtained from H{P, Q) is a constant of the motion and can be made 
equal to the initial energy. 

6. As the coefficients for friction or potentials are made to vanish, Q, Q, P, L(q, q, t), L(Q, Q), 
H(p,q,t) and H{P 1 Q) collapse smoothly to the expected simpler forms. 

7. Q, Q as well as q, q are related through simple differentiation with respect to time: Q — ^2; 
a = & 

9 dt ■ 



Previously derived examples for H exp include the damped free particle and the damped 
harmonic oscillator pQ. For a damped oscillator, with q + jq+ Lo 2 q = 0, the Q, L's, H's are: 

Q = ?exp(^) 

L(Q,Q) = f (Q 2 - (uj 2 - t)Q 2 ) 
H(P,Q) = f(^P 2 + (^-^)Qi) 

L(q,q,t) = f(q 2 +- 1 qq+^-q 2 ~u J 2 q 2 )exp( 1 t) 
H(p,q,t) - f(Jpp2 + l.pq + uj 2 q 2 )exp{yt) 

Solving the associated equation of motion for p(t) and q(t) and inserting the results back into 
H(p, q, t), shows that H(p, q, t) is a constant of the motion. H (p, q, t) is also equal to the initial 
energy under the right inital conditions, such as go = or go = for example. 

Q's, L's and H 's that satisfy conditions 1-7 for a particle moving in a constant gravity field 
or against friction for the lowest powers of velocity are shown in tables 1-a and 1-b. 

Table 1-a. 



eqn of motion 


Q definition 


L(Q,Q) 


H(P,Q) 


q + g = 
q + V = Q 

Q + 7<7 = 
q + bq 2 = 


Q = q 

Q = gexp(0) 
Q = gexpQi) 

Q = i(l - exp(-bq)) exp(bq) 


f(Q 2 ~2gQ) 
f (Q 2 -2vQ) 

f(Q 2 + ^Q 2 ) 

f Q 2 


f(-^P 2 + 2gQ) 

f(^P 2 + 2 M Q) 

ml \ -pi i 2 rfl\ 
2U 2f 4 > 
1 p2 

2m 



$8L represents the variation of the action integral S f 2 Ldt and is calculated from the 

Euler-Lagrange relation ^ ~ §§ = 0. 

''The need for condition 4 is seen for instance, whenever viscous friction is present. For 
quadratic friction, where time does not appear explicitly, H can indeed be found by using the 
Legendre transform in g-space. 
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Table 1-b. 



eqn of motion 


L(q,q,t) 




q + g = 
q + fi = 

q + iq = 
q + bq 2 = 


fif-Zgq) 

f (g 2 -2 W )exp(0) 

t(<? 2 + iqq + "2"9 2 ) exp(7i) 
fg 2 exp(26g) 


f (-W + 2 5 g) 

f (J I p 2 + 2 W )exp(0) 

f(^ 2 + ^w)exp( 7 t) 
exp(26g) 



The entries of table 1 have additional noteworthy characteristics. Each g-space L and H 
contains a dimensionless expansion factor 11 as indicated by Yan 0: 

X = exp[- J ^dt] (2) 

where G is the right hand part of the generic differential equation, q — G(q 7 q,t). Conservative 
potentials that exist in physical space are required by condition 6 to be explicity expressed 
within L(Q, Q) and H(P, Q) as well. Some cases, like the viscous example, require potential-like 
terms such as ^-Q 2 in L(Q, Q) and H(P, Q) for the conditions of H exp to be fulfilled. However, 
for the quadratically damped case, no such terms are needed. As will be seen, other frictional 
systems also require "pseudopotentials" that have the form f(Q, £oj £l> • • • j 4>0i 01) • ■ •) where Q 
is the transformed position and & and 4>i are the relevant coefficients of friction and potential. 

The tables 2-a and 2-b list the Q's, L's and H 's derived for various combinations of friction 
and potential. The symbol a is used to represent gravity g or dry friction [i that opposes the 
motion. Note that each L(Q, Q) and H(P, Q) contains some form of pseudopotental. 



Table 2-a. 



eqn of motion 


Q definition 


£(Q,Q) 


(P, Q) 


q + 7<j + a = 

q + bq 2 + 7<j = 

q + bq 2 + a = 
q + bq 2 + LU 2 q = 


Q = (g+M_ 2 « )eX p(f) 

-2^cosh(f ) + 4^ 

Q = i(l-exp(-6q)) 

x exp(6<7 + 2f) 

Q = \ (! - cx p(- & <?)) cx p(M 

Q = i(l - exp(-bg)) exp(fcg) 


f(Q 2 + ^Q 2 -2«Q) 

f(Q 2 + ^Q 2 ) 

f (Q 2 - abQ 2 - 2aQ) 
f(Q 2 

-^ln(&Q + l)(6Q + l) 2 

+^(&g + i) 2 -^) 


f (^P 2 -^Q 2 + 2aQ) 
f(^P 2 -^Q 2 ) 

f (^P 2 + a6Q 2 + 2aQ) 

m / 1 p2 
2 ff r 

+ £rhi(&Q + l)(6Q + l) 2 

-^(&g + i) 2 + ^) 



^Also known as an integrating factor. 
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Table 2-b. 



eqn 


of motion 


L(q,q,t) 




9 + 


79 4- 


-a = 


f((q 2 +iqq+^q 2 

-2aq- ^q + 2^ -2^-t + aqt 

+ajqt + sJ_) exp(7t) 
+-9 ^ ^— + -2 cosh(7<)) 


- -S-p - 2^ + sJ. + %pt) exp( 7 i) 

111 f I £ III ' » ' ' 

+ag+^ + 2^ + s_*) 


9 + 


bq 2 - 


f 79 = 


f (g 2 exp(26<? + 7 i) 
+i?(exp(6(?) - 1) cxp(feq + yt) 
+^(exp(6g) - l) 2 cxp( 7 t)) 


f (^ 2 exp(26g + 7 <) 

+^bP ( ex P(M - !) exp(6g + 7t)) 


q + 


bq 2 - 


fa = 


f ((g 2 - s)exp(26g) + f) 


f((^P 2 + f)exp(26 9 )-f) 


9 + 


bq 2 - 


f cj 2 (7 = 


f ((9 2 -^9+^)exp(26g)-^) 


f((^ 2 + ^-^)cxp(26g) + ^) 



It is readily demonstrated that both 5L(q,q,t) — and SL(Q,Q) = produce the correct 
equations of motion. Solving for g(i) and g(t) and inserting the results back into ff(p, 9, i), also 
shows that the H(p,q,t)'s are constants of the motion. Under requisite initial conditions, H is 
equal to the initial energy as well. Other combinations are also possible. 

There are general characteristics of solutions that are observed. The combination of two or 
more L(q,q,t)'s with differing friction terms results in their respective expansion factors being 
multiplied together pjj EJ ■ The Q _s P ace L f° r a combination includes the potentials and 
pseudopotentials of the simpler Us. For combinations of \x or g with a form that has viscous 
friction, the Q-definition contains an additonal fit) that is to be determined. 

As has been mentioned, pseudopotential functions in Q-space are seen to be important in 
finding solutions to dissiptive L's and -ff's in the Schuch formalism, but there is not yet a proved 
method for their determination. The current examples follow the provisional rule that if such 
a function can be found which satisfies conditions 1-7, then it is included. In some cases, only 
one simple possiblity exists that has the correct units, i.e. energy. For the free quadradically 
damped particle there evidently are none. In the case of the quadratically damped oscillator, a 
pseudopotential may be determined by integrating the equation of motion Q-space. Examples are 
presented to illustrate in more detail how the Q definitions and subsequent L and H combinations 
were constructed. 

Example 1: q + jq + a = 0. The ansatz Q-definition is taken to be: 

Q = «exp(2t) + /(*) (3) 

Since the invariance requirement of condition 2 means that the L must produce the correct 
equation of motion regardless of whether it is in Q or q space, Q-space may be used, which is 
more convenient. An initial L(Q, Q) is formed from the terms that correspond to each individual 
frictional coefficient and potential as listed in table la: 

2 

L(Q,Q) = j(Q 2 + ^Q 2 ~2aQ) (4) 

5L(Q, Q) = is calculated: 

SL(Q,Q) = -(2Q-^Q + 2a) = (5) 
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The Q above is differentiated to find Q and Q, which are used to convert the differential equation 
to g-space: The result is: 

2 

(g + 7 g)exp(|t) + /-^/ + a = (6) 
The following differential equation in / must therefore be solved: 

2 

/-!./ + a = aexp(|«) (7) 
The general solution for / is of the form, 



T —■~y q/£ '•y q; 

/(*) = a expUt) + c 2 exp(— h) + — exp(^i) + 4— (8) 
2 2 7 2 7 Z 



Using series expansions for c\ exp(^i) and c 2 exp(-^-t), then applying condition 6, it is found to 
be necessary that: 

ci = -3-5S c 2 = 5- (9) 



A valid Q definition is therefore: 



„ a at. .7 , a .—7 , a 
Q = (q - 3— + — expUt) - — exp — k) + 4— 10 
7 Z 7 2 7 Z 2 7 Z 



After combining terms, Q with resulting L's and H's are as shown in tables 2-a and 2-b. 

Example 2: q + bq 2 = 0. Derivation of the Q definition can precede if one notices that 
-^.(qexp(bq)) produces the desired equation of motion. It is possible to set up a differential 
equation for an unknown function / based on the assumptions: 

Q = fexp(bq); Q = f exp(bq) + bfqexp(bq) =qexp(bq) (11) 

The resulting differential equation is: 

/ + (&/- l)g = (12) 
The solution for / in terms of q is then: 

/=i(l-exp(-6g)) (13) 

A valid Q definition is therefore: 

Q= i(l-exp(-69))exp(6g) (14) 

The method above works for quadratic friction, but there is evidently no current technique for 
determining Q definitions of orders q 3 and higher. 

The L(q, q, t) that is obtained for Example 3 and shown in table 1-b is essentially the same as 
what is found in Havas Tartaglia pjj and Razavy |3]- Since time does not appear explicitly, 
the L's still match even if potentials are included [2]]. However, H exp requires by conditions 3 
and 5 that a system in real space be based on the kinetic momentum as opposed to the conjugate 
momentum. Because the resulting H(p,q,t) for this example uses p = mq instead of p — it 
differs from the i/'s of the aforementioned authors. 
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Example 3: q + bq 2 +-fq = 0. The initial Q-definition is assumed to be the Q for the highest 
order of friction is multipied by the lower order expansion factor: 

Q = 1(1 - exp(-bq)) cxp(bq) exp(^i) = ^(cxp(bq) - 1) exp(Ji) (15) 

A combination L(Q, Q) that includes the terms from each corresponding simpler L is assumed: 

HQ,Q) = j(Q 2 + ^-Q 2 ) (16) 
Then SL(Q, Q) — is calculated: 

5L(Q,Q) = ^ (2Q- yQ) = (17) 

Converting the differential equation to g-space obtains: 

(q + bq 2 + 7 <z) exp(6g) cxp(|t) = (18) 

The expansion factor is divided out and the desired result is confirmed with no other steps 
needed. 

Example 4: q + bq 2 + uj 2 q = The initial Q-definition is assumed to be: 

Q=^(l-exp(-bq))exp(bq) (19) 

The ansatz L(Q, Q) should include uj 2 Q 2 but it will be more convenient to leave this term out 
for the time being: 

L(Q,Q) = ^-(Q 2 ) (20) 

The 5L(Q,Q) is calculated: 

777 

SL(Q,Q) = -(2Q) = (21) 

Converting the differential equation to g-space: 

(q + bq 2 ) cxp(bq) = (22) 

It is clear that an u 2 qexp(bq) term is required to complete the equation of motion. Even 
though friction is present, there is no explicit time dependence. It is possible to work backward, 
converting the needed term to a function of Q by using the Q-definition. 

LO 2 

cu 2 qexp(bq) = —\n(bQ + l){bQ + l) (23) 
o 

The right hand side of equation (23) above is integrated to complete the desired L(Q, Q) which 
results in the complicated looking expression at the bottom of the third column in Table 2a. 
By expanding these terms into a power series it is seen that the conservative potential oj 2 Q 2 is 
present. A constant —552- has also been added to satisfy condition 6. For q + bq 2 + a = the 
same procedure as just described in this example can be used. 

Conclusion 

It is possible to construct Lagrangians and Hamiltonians with additional combinations of 
friction and potential terms that are consistent with a previously published Lagrange-Hamilton 
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formalism for dissipative systems. The resultant L's and Lf's offer new possibilities for canonical 
quantizations of frictional systems. Other areas of future study include the determination of a 
general method for finding Q definitions for higher order frictional cases. The physical signif- 
icance of the L's and H's, such as may exist, remains to be completely explained. A proof of 
uniqueness for L's and LT's may require additional, as yet undefined conditions to ensure unique 
physical interpretations. 
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